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Abstract 

Relativistic heat transport in electron-two-temperature plasmas with density gradients has been 
investigated. The Legendre expansion analysis of relativistically modified kinetic equations shows 
that strong inhibition of heat flux appears in relativistic temperature regimes, suppressing the 
classical Spitzer-Harm conduction. The Seebeck coefficient, the Wiedemann-Franz law, and the 
thermoelectric figure of merit are derived in the relativistic regimes. 
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The relativistic effects of hot electrons in laboratory plasmas have attracted much inter- 
ests in the past few decades 111, particularly, in the context of the current drive mode in 

nn 

tokamaks and the electron cyclotron heating in various confinement devices |2J, |3(. In the 
theoretical arena, Braams and Karney (BK) first presented the relativistic Fokker-Planck 
equations with the extended Rosenbluth potentials |4( , and applied the equations to deriva- 
tion of the relativistic electrical conductivity 5|. Using the BK collision integrals, Shoucri 
and Shkarofsky developed a numerical code to survey the relativistic effects on electron 
cyclotron wave, fast wave, and lower hybrid current drive mode in tokamaks [6]. Making 
use of the Chapman- Enskog expansion, Mohanty and Baral derived relativistic transport 
coefficients including magnetic field effects 

Rapid heating of plasma often leads to a bi-Maxwellian electron distribution, consisting of 
a bulk and high-energy tail, while maintaining relatively cold ions. In fact, the appearances 
of two-temperature spectrum for electrons have been observed in some experiments with 
high-intensity lasers |<a]. The well-pronounced tail and its velocity moment, which determine 
the transport properties, are in the relativistic temperature regime and, therefore, important 
to the fundamental study of relativistic electron transport. 

In this Brief Communication, the relativistic transport theory presented in the previous 



paper 



is expanded, aiming at the numerical simulation of high-temperature ignition 
plasmas, and celestial plasmas. We focus on a problem relevant to heat flux inhibition 
due to relativistic effects of electrons within the framework of the relativistically corrected 
Spitzer-Harm (SH) formula for electron-two-temperature plasmas with density gradients. 
The formula is fully consistent with the current-neutral condition, so that one can readily 
couple the transport coefficients with fluid codes [lol |. 

The relativistic thermal conductivity is derived below, along the manner developed by BK 
. Begin with the Legendre expansion for electron distribution function, viz., f(r, p, t) ~ 
fo(r,p,t) + (p/p) • fi(r,p, t) + hiqher order terms for the small parameter which is related 
to a characteristic field strength 9]. Introducing the relation of p x = pcos<f), a relativistically 
extended kinetic equation is averaged over the solid angle Q, i.e., < • • • >= J ■ ■ -d(fi/47r). 
After the manipulations, we obtain the first order equation in the form of 

dfi x _e^dfo_ = f5fiA 



dt dx m c dfi \ 5t 
where v = c/i/T, /i = p/m c, and T = y/l + /i 2 . Using the transfer cross section of the rela- 



tivistic Mott scattering a l ei1 the collision term of Eq. (JH) can be approximated by (Sfi x /St) c ~ 
—niV<j\J\ x = —feifix, where n.j is the number density of ions, v ei = (niY ei /c 3 )(T/ fi 3 ) is the 
electron-ion collision frequency, Y ei = 47r[Ze 2 /(47reomo)] 2 lnA, and m , Z, and InA are the 
electron rest mass, the averaged charge number, and the Coulomb logarithm, respectively 

a. 

In Eq. the effects of magnetic fields are ignored. This approximation is valid for 
v e i ^> u c , where uj c = \eB\/(Tmo) is the electron cyclotron frequency. The validity condition 
gives the allowable parameter range of the magnetic field strength of 

|B| « 1.7 x 10' r * ( f "■ 3 ) (—) G. (2) 
I I (p2 _ i) 3 / 2 V10 27 cm- 3 / V !0 / 

In highly compressed targets irradiating by a relativistic laser pulse the dense plasma 
parameters are typically (T — 1) ~ 10 _1 , rij ~ 10 26 cm -3 , Z rs 3.5 (carbonized deuterium- 
tritium), and InA « 5. For such parameters, the right-hand side (RHS) of Eq. (j2J) reads 
about 130 MG. Around the tenuous coronae, the laser pulse drives relativistic currents, and 
induces the self-magnetic fields of magnitude B ~ m Q u) pe c/e ~ 10 2 MG. Recent numerical 
simulation indicates that at the surface, an intense magnetic field of B < 280 MG prevents 
hot electrons from penetrating into the higher density region Q|. The electrons near the 
channel envelope where the magnetic field is strongest, as well as lower energy electrons, 
tend to be magnetically trapped [13], since their Larmor radii are comparable to or even 
less than the channel radius. This stopping effect seems to be subject to the Alfven current 
limit which is irrelevant to the limit of energy flux. For a self-focusing electron beam, 
a fraction of energetic electrons, running along the channel axis, cannot be trapped 
and generate the relatively small magnetic fields of B < 10 2 MG in the denser plasma, as 
was shown in Ref. The penetrating electrons have a largely anisotropic momentum 

distribution, which cannot be treated by the diffusion approximation employed here. In a 
highly compressed region, however, the transported electrons are expected to be thermalized 



via dissipative processes 



14, 



, 3|, and the beam type transport may become diffusive, further 
decaying the magnetic fields. Although, diffusive transport plays a significant role in heating 
the final compressed fuel, the details are still not well understood. Hence, here we investigate 
the fundamental transport properties in the parameter regions of the highly compressed 
ignitor plasma, where magnetic field effects can be fairly neglected, as far as Eq. (J2J is 
fulfilled. I also mention that the density gradient of ablative plasma is likely to be steep in 



the higher density regions 10], so that nonuniform effects are taken into account here. 

For a quasisteady condition of df\ x /dt ~ in Eq. (JTJ), i.e., omitting the electron inertia, 
the anisotropic component of electron distribution function is given by 

f» M *'(fM^-?°*>). (3) 
riiY ei y 1 + y 2 dx m c 2 y/l + fj, 2 o/j, J 

Heat flux of relativistic electrons can be defined by q x = niQC 2 J °° j (T — l)v x fi 2 fdQdfi. 
Integrating over solid angle, yields q x = |7rm c 3 / °° fi x fi 3 (T — l)/Tdfj,. Making use of Eq. Q, 
this may be written as 

4nm c 7 ( f 00 J 0/o d/i eE x ^ /i 6 0/ o ^ 



3niY ei \J Q 1 + /j 2 dx m c 2 J ^1 + jjfi dfx 

The longitudinal electric field E x in Eq. can be determined by the current-neutral con- 
dition j x = — e / °° J v x fi 2 fdVldfi = —^nec J °° fi x fi 3 /Tdfi ~ 0. That is, 

eE x _ Iq (i + ^)3/2 d/i. 

For the isotropic component f (x,fi), I employ the superposition of the two-temperature 
populations of electrons, 

r , \ 1 tt-^ n e Ax)aAx) r , . r -~| ,„ s 

MW) = K2[a . {x)] exp [- % (x)^lT^j , (6) 

where K u (aj) is the modified Bessel function of index v with its argument of <x,(x) = 
ttiqc 2 /Tj(x), and j — c,h indicate the cold and hot components, respectively. The normal- 
ization is given by n e (x) = Zrii(x) = J2j n e,j( x ) — 4vr fo(x, /i)/x 2 d/i. 

In a steep temperature gradient plasma, depending on the collisional mean-free path, 
A, the transport properties may not be locally defined. In this sense, local transport 
theory is valid only for the case of A \Lt\, where Lt = T/(dT/dx) is the charac- 
teristic length of the temperature gradient. Concerning the relation of \dT/dx\ ~ e\E x \ 
derived from Eq. (jSJ), the parameter range involving the electric field can be estimated as 
\E X \ < E c ~ 10 12 (100 keV /T)(Z 2 rii/l0 27 cm- 3 )(lnA/10) V/m. In the case of ignitor physics, 
the relativistic electron transport establishes the temperature gradient in the high-density 
plasma. Assuming the spatial gradient of AT/Ax ~ —100 keV/100 /im, the electric field 
strength can be estimated as \E X \ ~ 10 9 V/m [see also Eq. (fTTjl below]. For T < 10 2 keV, 
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Z 2 m ~ 10 27 cm" 3 , and InA « 5, we read \E X \ < 10~ 2 £ c . For the case of \E X \ > (0.01-0. l)E c , 
one may solve kinetic transport equations to determine the full self-consistent spectral dis- 
tribution, instead of using Eq. (jHJ) [16]. 

Substituting Eq. (JHJ) into Eqs. (@J) and (0), we obtain the relativistic heat flux of 
q x = —K re [(dT h /dx) for the temperature gradient of hot electrons, and may decompose 
the coefficient as n re i = f re i ksr. Here, K SU (T h ) = 256(2ir) 1/2 e 2 ) T^ /2 /(Ze 4 ml /2 \nA) is the 



familiar nonrelativistic SH heat conductivity of the Lorentz plasmas 



17j , and the factor f re f 



corresponds to the relativistically corrected flux limiter which can be expressed as 

frel = ^^«I /2 {Cl,e + OC^h + € [C^Q^a.) + C^K)]}, (7) 

where the abbreviations are 

0C 2 , C [Qi(a c ) + C 3 , c e 2 (a c )] + C %h [@i(a h ) + C 3 , h e 2 (a h )] 



a c C 2 , c 02(a c ) + a h C 2 hQ 2 (a h ) 



(8) 



c . = C2 J -[C73 J -ei(a J -)-e 3 (a J -)] (9a) 



.1 



Ou = ^^H, (9b) 

n e K 2 {a j ) 



K 2 (a j ) 



^ , n ( , 1 2 42 120 120\ 

©i(«i) =1 - — + 72 + 73 + TT + ^ exp(-a,) + % Ei(-a 

y CXj j j j j J 

( 1 2 6 24 24 \ 
2 (a i ) = ( 1 - — + -f--3--4--5 ex P(-«i) + %Ei(-o; j 

Y Uj Ctj Ctj Ui'j U'j J 



(10a) 



(10b) 



^ . , /48 288 720 720\ , . . 

©3(«i) = (^2+^3- + 7X + 7X ex P(-";)> ( 10c ) 

where Ei(— atj) is the exponential integral function, and 8j = d\nn e ^/d\nTj and 6 = 
d\nT c / dlnTh reflect the nonuniformity of plasma. Namely, for 6^0, the formula describes 
the energetic transport in the plasma that the cold electron component is isothermal, and 
for Sj — + and —1, in the plasma that the electron component j is isochoric (n e j = const) 
and isobaric (n e jTj = const), respectively. 



The geometrical constraint of Vn e j \\ VI} due to ignoring two-dimensional (2D) effects 
means that thermoelectric magnetic fields, which can be prominent, for example, in intense 
laser-plasma interactions [3], are not taken into account at the moment. 2D effects are 
important, because they prefer to short out electric fields and pinch directional flows by 
the toroidal magnetic fields. The complexities of magnetic inhibition in heat flux might be 
effectively considered by introducing a reduction factor /b < 1: Bohm's /b = (1 + ^ c t) _1 
or Braginskii's /b = (1 + u; 2 t 2 ) _1 18|, where r denotes a collision period. That is, one 
can practically utilize the cross- field conductivity approximated by n±_ re i ~ /b frei k sh- Note 
that Eq. (J2J) reflects the much smaller Hall parameter uj c t <C 1, such that /b — > 1, and 

^_I_reZ — ^||reZ ^rel- 

In the following, more elemental issues are investigated, i.e., relativistically extended lon- 
gitudinal thermoelectric effects. With regard to the longitudinal thermal diffusion that devel- 
ops an electrostatic potential, one should note the important relation e = LT h [eE x / '(moc 2 )]. 
For the special case of n e ^/n e — > 1, C2 jC — > 0, and 9 —>■ 1, namely, the one-temperature 
model for electrons, the self-consistent electric field Eq. (JHJ) reduces to 

I 



e ~ — 

a 



— — - + C 3 [a,d) 
B 2 a 



(11) 



where a = a.,- and C3 = C-^j. In the thermoelectric point of view, the relativistic Seebeck 
coefficient can be defined by s = ae/e. The temperature dependence of Eq. (fTT|) is shown 
in Table U for 5 = 5j = and —1. In the nonrelativistic limit of a ^> 1, Eq. (1111) for 5 = 

n 

asymptotically approaches e — > — 5/(2a) [9|. In the isobaric case of 5 = —1, owing to the 
pressure-balance effects, the field strength reduces to 50 — 60 % of the 5 = case. Noted is 
that in this case the flux limiter j re \ does not depend on 5, and the similar property appears 
again in the following other cases. 

In Fig. the temperature dependence of the relativistically corrected flux limiter are 
shown. The ratios of hot/total electron density are chosen for n e ^/n e = 0.1 — 1, fixing 
the temperature scale length equal, 6=1. For the electron- two-temperature models of 
n e ,h/ n e 7^ 1 (C 2iC 7^ 0), set the temperature of cold component to a c = 10 2 (T c = 5.11 keV) 
as an example. The densities are set to be uniform (5j = 0), except for the case of n e ^/n e = 
0.1 that the nonuniformity (5j = —1) is taken into consideration. Actually our major 
interests are in the relativistic heat flux carried by the high-energy tail electrons of T > 
To = ( a h — Q ; c)~ 1 l n (C2,c/C , 2,h), where the spectral population of hot electrons is larger than 



that of cold ones. As expected, in the lower energy regions of V < Vq, energy transport by 
cold components is dominant. In Fig. [TJ such criterion seems to appear as pseudo cut-off in 
the lower temperature region. 

Now one finds that the heat flux is strongly inhibited in the relativistic regime. For 
example, in the one-temperature model for electrons, the flux limiters are f re i ~ 0.73 for 
T = 0.1m c 2 (a = 10) and f re i ~ 0.37 for T = m c 2 (a = 1), as shown in Fig. [T] (solid 
curve). This is due to the drift velocity carrying heat asymptotically close to the speed 
of light. Moreover, it is found that a fall in the hot electron population leads to further 
decrease of the conductivity, and indeed, the degree of the depletion reflects the abundance 
of hot electrons. Regarding the electron transport in laser-produced plasmas, typically a 



consistent with 



flux limiter of order of 10 -2 — 10 _1 has been empirically employed 

the experimental results [3]. In this aspect, the present results imply that the relativistic 
effects on sparsely populated high-energy tails can also participate in lowering the flux limit. 
These properties do not largely depend on 5j as seen in Fig. [T] For example, in the case of 
n>e,h/n e = 0.1, the difference of the flux limiter between the case of 5j = (dotted curve) 
and —1 (crosses) is about 10% at most. 

Taking the limit of n h ^ e /n e — > 1, C 2<c — > 0, and 9 — > 1, Eqs. (fTjl- lfTUj) reduce to 



rel 



(27T) 1 / 2 a 7 ' 2 



e?(q) 
e 2 (a) 



+ e 3 (a) 



(12) 



384 K 2 (a) 

This corresponds to the standard relativistic SH heat conductivity having the temperature 
dependence of khm{T) = f re i(T)KsE.(T) oc T 2 — T 5 / 2 (<]], which exhibits the asymptotic 
property of ksh(^) oc T 5 / 2 (f re i — > 1) in the nonrelativistic limit of a > 1 whereas 
KDT (T) = [5(27r) 1 / 2 /32]a 1 / 2 Ks H (T) cx T 2 (f rel oc T" 1 / 2 ) by Dzhavakhishvili and Tsintsadze 
in the ultrarelativistic limit of a -C 1 |20j]. These characteristics are also shown in Fig. ^ 
(solid curve), and summarized in Table [I] 

Here let us take the ratio of the thermal to electrical conductivity. The key relation 



is known as the Wiedemann- Franz law for metallic states of matters |2l| . The ubiquitous 
nature is derived from a simple assumption of the elastic scattering of conduction electrons. 
As for fully ionized plasmas, the relativistically extended law can be expressed as 

«hm T Q?(a) + 6 2 (Q)9 3 (a) ^ 
CTbk e 2 Q 2 2 (a) 

for the case of n e)h /n e = 1. Here, a BK {a) = -[{2it) 1 ' 2 /96][a 7 / 2 6 2 (a)/ K 2 {a)]a s > 0, 



and cs(T) = 64(27r) 1 / 2 e 2 1 T 3//2 / (Ze 4 7?\/ 2 lnA) stands for the nonrelativistic Spitzer conductiv- 
ity. Evidently, the ratio depends on the temperature only, without involving the intrinsic 
parameters of plasmas. Figure El shows the temperature dependence of Eq. (jl3j) . In the 
nonrelativistic limit of a 1, it asymptotically approaches khm/o"bk — ► ^sh/os = 4T/e 2 
Q. This value slightly decreases as the temperature increases, to take the minimum value 
of (ftHM/cBiOmin = 3.92T/e 2 at a = 19.6 (T = 26.1 keV). As seen in the figure, it increases 
up to kdt/cbk = 5T/e 2 in the ultrarelativistic regime. It may be instructive to men- 
tion that the transport equation of the Fermi liquid in metals or condensed plasmas yields 
k/o ~ 7r 2 T/(3e 2 ) ~ 3.3T/e 2 J2]} 22^ which is lower than (/tHM/omOmin in the ordinary 
plasmas. 

The heat conductivity holds the larger power index of temperature. Thus, fast heating 
of plasma can drive the nonlinear heat-wave accompanied with a well-defined wave front, 
where an electrostatic field tends to be well developed. This leads to an idea that such a 
thermally non-equilibrated plasma can be essentially compared to a thermoelectric converter. 
And, in general, its efficiency can be quantitatively evaluated by invoking a thermoelectric 
figure of merit. Along the conventional notation used in material physics, we now define 
the thermoelectric figure of merit by Z = s 2 ctbk/«hm = {cue/ 'e) 2 (<7BK/ 'khm) for n e ^/n e = 1. 
Making use of Eqs. (jllj) and ([13)1 . this multiplied by T can be written in the dimensionless 
form, 

[Q 1 (a) + e 2 (a)C 3 (a,5)} 2 

e 2 (a) + e 2 («)e 3 («) 1 ) 



Note that Eq. (jT4|) depends on 5, in contrast to Eqs. ([12]) and (jl3j) . For 6 = and — 1, 
the temperature dependence of Eq. (|14|) and the coefficient s are also shown in Fig. |21 
It is found that for a 1, Eq. (jT3|) asymptotically approaches ZT(5 = 0) — > || and 
ZT(5 = -1) -> whereas for a < 1, ZT(5 = 0) -> f and ZT(5 = -1) -> |. Particularly, 
in the nonrelativistic plasmas with uniform density, i.e., a ^> 1 and 5 = 0, one can extract 
the higher figure of merit ZT = 1.56, compared with the typical thermoelectric materials 
which provide ZT ~ 0.4 — 1.3 as indicated by the arrow in Fig.lJlj^]. Notice that the Carnot 
efficiency can be achieved for ZT ^> 1. The dimensionless values of (KHM/c"BK)(e 2 /T), (se) 2 , 
and ZT for some a and 5 values are summarized in Table |HI 

In conclusion, I have derived solutions for the heat conductivity and related thermoelectric 
coefficients in a relativistic nonuniform plasma. These results indicate that the relativistic 



S 



effects on the high-energy tail electrons significantly limit the heat flux. This mechanism 



might p 
physics 



ay an additional role of the stopping of relativistic electrons in the context of ignitor 
15j . although this work ignores 2D thermoelectric effects such as V±n x VnT, which 



may be important for typical ignitor geometries. 
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Hot electron temperature T h lm Q c 2 



FIG. 1: Hot electron temperature dependence of the relativistic Spitzer-Harm (SH) heat conduc- 
tivity normalized by the nonrelativistic one: The flux limiter of hot electrons f re i(Th) = n re i/ '«sh 
is shown for density ratios of n ej /i/n e = 1 [Eq. (|12j) : solid curve], 0.5 (dashed curve), 0.2 (dotted- 
dashed curve), and 0.1 (dotted curve). For the case of n ei h/n e ^ 1, T c /niQC 2 = 0.01 and n e j = const 
are chosen as examples. For comparison, we plot another case of n e j oc TJ , only for n e /j/n e = 0.1 
and T c jmQ(? = 0.01 (crosses). 
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FIG. 2: Electron temperature dependence of the dimensionless thermoelectric figure of merit ZT 
(solid curves), the Wiedemann-Franz ratio khm/cbk multiplied by e 2 /T (dotted-dashed curve), 
and the Seebeck coefficient s multiplied by |e| (dotted curves). Note that in the present case 
tHivi/ffBK is independent on 5 and takes the minimum values at T/rriQC 2 ~ 0.051. 
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TABLE I: Temperature dependence of f re [ and ae for 5 = and —1. 

a T (keV) ae|<5=o ae|«5=-i frei a 

<1 -2-1 

1 5.11 x 10 2 - 2.0178 -1.0178 0.36792 

5 1.02 x 10 2 -2.1692 -1.1692 0.62668 

10 51.1 -2.2692 -1.2692 0.73598 

20 25.5 -2.3560 -1.3560 0.83044 

100 5.11 -2.4638 -1.4638 0.95529 

> 1 -2.5 -1.5 1 

a frei = khm/ksh Q to give f re i ~ 1 an( l ^ T~ 05 [^(J for q » 1 and <C 1, respectively. 
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TABLE II: Temperature dependence of (achm /f bk ) (e 2 /T) , (se) 2 , and ZT for S = and — 1. 



a: 


T (keV) 


(«HM/^BK)(e 2 /T) a 


{se) 2 \s=o 


(se) 2 | 5= _i 


ZT\s=o 


ZT\ &= _ 


< 1 


- 


5 


4 


1 


4/5 


1/5 


0.5 


1.02 x 10 3 


4.9032 


4.0145 


1.0072 


0.81874 


0.20542 


1 


5.11 x 10 2 


4.7331 


4.0717 


1.0360 


0.86026 


0.21888 


5 


1.02 x 10 2 


4.0886 


4.7053 


1.3670 


1.1508 


0.33434 


10 


51.1 


3.9529 


5.1495 


1.6110 


1.3027 


0.40754 


20 


25.5 


3.9221 


5.5509 


1.8388 


1.4153 


0.46883 


100 


5.11 


3.9669 


6.0704 


2.1428 


1.5303 


0.54170 


> 1 




4 


25/4 


9/4 


25/16 


9/16 



q cbk(T 1 ) is introduced in Ref. 
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